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We present detailed calculations for the Casimir force between a plane and a nanostructured 
surface at finite temperature in the framework of the scattering theory. We then study numerically 
the effect of finite temperature as a function of the grating parameters and the separation distance. 
We also infer non-trivial geometrical effects on the Casimir interaction via a comparison with the 
proximity force approximation. Finally, we compare our calculations with data from experiments 
performed with nanostructured surfaces. 
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1. Introduction 

The Casimir interaction at finite temperature consists 
of two parts: a purely quantum one which subsists as 
T — K involving the zero-point energy of the electro- 
magnetic vacuum IJ and a thermal one [2^ which takes 
into account the real thermal photons emitted by the 
bodies. For the thermal part to become noticeable, fre- 
quencies relevant to the Casimir interaction must fall in 
the range of relevant thermal frequencies. For this reason, 
at room temperature the thermal part of the Casimir in- 
teraction becomes important for separation distances of 
the order of microns or tens of microns. At those sepa- 
ration distances, the absolute value of the Casimir force 
is small as it decreases as the inverse third power of the 
separation distance. Experimentally, there is thus a trade 
off when one wants to measure the thermal component 
of the Casimir interaction: at small separation distance 
where the Casimir force is comparably large, the ther- 
mal part is small whereas at large separation distance 
where the thermal part is large, the total Casimir force 
is small. A solution is to use "large" interacting bodies 
to maximize the Casimir force at large distances 

The above reasoning is well adapted to the parallel- 
plates geometry (or a plate-sphere situation when the 
radius of the sphere is large) where the only character- 
istic length is the separation distance. In the case of a 
plate-grating situation additional characteristic lengths 
such as the grating period or the corrugation depth are 
to be taken into account. The importance played by the 
thermal part of the Casimir interaction is highly non- 
trivial in the case of the plate-grating geometry and full 
calculations become necessary. First calculations for this 
geometry for perfect reflectors at zero temperature have 
used a path integral approach [3]. A first exact solu- 
tion for the Casimir force between two periodic dielectric 
gratings was given in (5] . Alternatively, the Casimir force 
in such geometries can also be calculated using a modal 
approach jB]. In this paper we use the scattering ap- 
proach to Casimir forces [THSj to calculate the Casimir 
interaction between a plate and a grating at arbitrary 
temperature. We study the contribution of the thermal 



part as a function of the grating parameters and assess 
the validity of the proximity force approximation (PFA) . 
We finally compare our results with experimental data 
presented elsewhere [TOl [11] . The most important result 
is that in the grating geometry the thermal contribution 
to the Casimir interaction is overall enhanced and occurs 
at shorter separation distance, which opens interesting 
perspectives for new experiments. 

2. Theory: Thermal Casimir force between gratings 

We study the Casimir interaction within the scattering 
approach between a plate and a ID lamellar grating as 
depicted in figure [T] Above the grating z > 0, we have 
a homogeneous region labelled I characterized by a per- 
mittivity e^. Below the grating z < —a, a homogeneous 
region labeled III characterized by a permittivity . The 
plate is characterized by a permittivity ep for z > L. In 
the grating region —a < z < 0, the permittivity is a peri- 
odic function of x, e(x). The Casimir force per unit area 
Fp^g{L) between a plane and a grating separated by a dis- 
tance L is calculated in the scattering formalism taking 
into account the finite temperature T 

FpjL;T) = 2TTkBT 

Y.' fftr ((1 - M„r^ dLMn) dk^dky 

where the prime on the sum means that the term with 
n = is to be multiplied by a factor 1/2. This ex- 
pression takes into account the contribution of a ther- 
mal field of real photons whose mean number n(w) = 
^^hu/ksT _ 2j follows a Planck's law to the contribu- 
tion of the virtual photons emerging from the electro- 
magnetic vacuum. In the following we will use c = 1 for 
convenience and work with a generalized complex fre- 
quency n — u + having real and imaginary parts uj 
and ^ respectively. 

The function Al„ is evaluated at the Matsub- 
ara frequency ^„ — ^^"^^^ and reads Ai{i^n) — 
Rp(i^„)e~'^^Rg(z^„)e~'^^. It describes a full round-trip 
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of the field between the two scatterers, that is a reflec- 
tion on the plate via the operator Rp, the free propa- 
gation from the plate to the grating corresponding to 
the translational operator e~'^^ with the imaginary wave 

1/2 

vector K = —i]iz = {^i^n + + ^y) ' scatter- 
ing on the grating via the reflection operator Rg and 
a free propagation back to the plate. The vector k^: 
of dimension 2N + 1 collects the diffracted wavevectors 
= (fc^ - N^, ■■■ ,k^,--- ,k^ + N^) where N is the 
highest diffraction order retained in the calculation. The 
integration in eq. ([T]) is restricted to the first Brillouin 
zone i.e. — ^ < fcx < f ■ 

The plate's reflection operator Hp is diagonal and col- 
lects the appropriate Fresnel reflection coefficients Rp — 

diag(^,|^) with .p = {epa+^l + klf'- 
We calculate the grating reflection operator Rg in 
the framework of the Rigorous Coupled Wave Analysis 
(RCWA). We use a method of resolution inspired by the 
formalism presented e.g. in [12) which we will briefly 
outline below. 

The physical problem is time- and j/-invariant. A 
global dependence in e*'^'^"^""*^ can be factored out of 
all the fields. 




region III: 



FIG. 1. A ID lamellar grating and its associated reference 
frame. 



For an incident wave characterized by a wave vector 
k^^^ = {ap,ky,—^p^) a grating structure with period d 
will generate an infinite number of reflected waves with 
wave vectors k^"^ — (a„, /cy, 7!*'') and transmitted waves 

with wave vectors k|"^ = (q!„, ky, — 7!* ■*). — k.^ + n"^ 
and7i'/'^ = (Q/,c.2-a2_fc2)i/2_ 

In region I, the y-components of the fields are written 
as a Rayleigh expansion involving incident and reflected 



fields of order p and n respectively 

+00 



Ey{x,z)= Snpli'^e^'^''''' 
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whereas in region III, the Rayleigh expansion involves the 
transmitted fields 
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Whether in region I or III, the x-components of the 
fields are obtained from the j/-components thanks to the 
Maxwell's curl equations 
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In the grating region —a < z < owing to the peri- 
odicity along the x direction the fields as well as the per- 
mittivity e{x) and its reciprocal l/e{x) can be expanded 
in Fourier series. We have 
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With those notations, we are able to express the 
Maxwell's curl equations in a compact matrix form. Let 
F be a column vector collecting the Fourier components 
of the fields F = (63,, e^,, h^., hj^)^, we may write 
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In the above equation, a — diag(a„), 1 is the identity 
and e, resp. 3, are ToepHtz matrices whose structure 
is defined as having elements {e„,n > 0} on the first 
fine and elements {e„,n < 0} on the first column. Note 
that in accordance with ref. |13j . we have replaced the 
matrix e by 3 in the lower-left block of the matrix M. 
This constitutes an improvement with respect to ref. [S] , 
where this replacement had not been done. The matrix 
F has dimension 4 x 1 in units of 2N + 1. A particular 
column of F corresponds to a particular incident order 
p. In the following, bold quantities are matrices whose 
dimensions will be given in units of 27V + 1 if not trivial. 
As an example, in eq. (6b) a, e, 3 and 1 are matrices of 
dimension 1 
2x2. 

In [5] eqn. 



1 whereas Mi and M2 are of dimension 



(6a 



had been numerically solved. Here we 
follow a different path which has proven to lead to more 
stable numerical calculations. Because of the block anti- 



diagonal structure of matrix M, eq. (6a) can be recast as 
a Helmholtz-like equation for the electric fields provided 
that M is independent of z which is the case for the 
lamellar gratings we consider here 



= M1M2 



This equation is solved as 
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^-VmJ^zq- TffT^Qj-Q c+ and C~ are unknown coefficients 
to be determined. Let (p, A be respectively the eigen- 
vectors and eigenvalues of the matrix M*^"^^ such that 
M*^*^) = ^diag(A) 0^^. Writing explicitly the expres- 
sion for e±\/'vi(=)2^ include the matrix in the 
unknown coefficients and C~; furthermore we want 
to avoid exponentially growing solutions at z = —a. Fol- 
lowing the prescriptions in [12[ we finally arrive at 



^ ) (z) ^e^^C+ + 0e-^(^+°)C- 



[j^^ (z) = M5-VVAe^''C+ 



(8a) 
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where eq. (8b) has been obtained by injecting eq. (8a 



into (^jj^j = M^^d, (^''j from eq. (6b I. 

We can use eqs. ([2| and eqs. ^ to write the fields at 
z = and eqs. ([3| to write the fields at z = —a. In 
compact matrix form, this leads to 
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and 



F(0) 



kyOL 




EiOj'^—k,^ 

1 




kyCX. 


EiUj'^—k'^ 




eiUj'^—k'^ 

1 


kyOt 




1 


^iUJ^-kl 






kyOt 


EiLO^ — k^ 




Cihj'^ — k^ 
1 



i/ie hih 



jfT— e 
■ja—h 



(10) 



JCT— e 
ja—h 



where 7^'/*) = diag(7i*^*'') and we have introduced the 
basis of polarizations a we use, denoted e and h. The 
polarizations a = e,h are defined by imposing Hy — 
and Ey = respectively. Hence, in the above equation for 
incident a = e waves, we impose 1°'='^ ^ 1 and 1'^='* = 
and vice- versa for incident a — h waves. Note that tg, 
t/,,, Tg and r/i are of dimension 2x2 whereas iee, ieh, ihe 
and ihh are of dimension 2x1. Other dimensions are 
deduced so as to be consistent with those of F. 

Evaluating eqs. dsl at 2; = — a and z = and iden- 
tifying with eqs. (pT and (10 1 leads to a linear system 
of equations of dimension 8{2N + 1) for the 8(27V -I- 1) 
unknowns C+, C", Rf'^', RC*), T^'^) and tC'). Neverthe- 
less, it is numerically more stable to eliminate the reflec- 
tion and transmission unknowns from this system and to 
solve instead a reduced system of dimension 4(2iV -I- 1) 
for solely C+ and C^. All done, this system reads: 



4 



0- 




(0 + r,r^^V 



jcr— e 
jcr— 



(11) 



where we have defined V = M^^0-\/A. Once the un- 
known coefficients C"*" and C~ are determined by solving 
eq. (11), the reflection and transmission coefficients are: 
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Resolution of eq. (11) and eqs. (12) ffist for incident 
a = e waves and then for incident a ~ h waves leads 
to the complete reflection and transmission matrices Rg 
and T„: 
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The force at T — K is recovered by the substitution: 
27rfcsT V' ^ /i / dC (14) 



We define two quantities, dpiyL) and r]p{L) to assess 
respectively the effect of the finite temperature and the 
deviation from PFA: 



Fp^g{L;T = 300 K) 
FpjL;T = OK) 



(15) 
(16) 



where F^g"^ 



{L;T) = i/;Fp,p(L(x);T)dx and 
Fp^p{L;T) the force between two plane-parallel plates is 
given by the Lifshitz formula |14| . 



3. Numerical evaluations 

We now study the thermal Casimir interaction between 
a gold plate and a doped silicon grating. Therefore, 

= esi{uj) and = eAui^)- The plate and the grating 
are separated by vacuum = eq = 1 Vw. As described 
in [llj , the permittivity of gold is taken from experimen- 
tal data extrapolated to low frequencies by the Drude 
model. The permittivity of doped silicon is modeled by 
a two-oscillator model: one describing the intrinsic part 



of silicon and the other one describing its metallic behav- 
ior at low frequencies [TS]. The metallic part of doped 
silicon is determined by a doping level of 2 x 10^^ cm""^. 
The situation is characterized by four length scales: the 
separation distance L, the corrugation height a, the grat- 
ing period d and the corrugation width w (for this last 
quantity, we will prefer to work with the filling factor 
/ = w/d). A complete analysis would in principle involve 
full Casimir force calculations in a four dimensional pa- 
rameters space. Instead we explore here the parameter 
space at fixed filling factor and grating period d corre- 
sponding to the experimental set-up in [llj . 

In figure[2]we illustrate the effect of the temperature in 




FIG. 2. Effect of temperature on the Casimir force given by 
the ratio -dp as a function of the separation distance L and 
corrugation depth a. 

the plate-grating geometry by plotting the ratio •dp{L^ a) 
as a function of both the separation distance L and the 
trench depth a. The grating period d and the filling fac- 
tor / are fixed respectively at d = 400 nm and / = 0.5. 
For this choice of materials, we find over the whole range 
of parameters 'i?i?(L, a) > 1 so that the thermal photons 
always lead to an increase in the Casimir force. For a = 
we recover the two-plate configuration and we have nec- 
essarily hmi_j.o t?_F(L, a = 0) = 1. The total temperature 
effect t?_F(L, a = oo) — i?F(i, a = 0) increases with larger 
separation distances L. The limiting value 'dpiL, a = oo) 
is reached for larger a as the separation distance L in- 
creases since this limiting value rather means a » L. 
Interestingly, for a fixed separation distance L, there 
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is a steep increase of -dp as a function of the corrugation 
depth a towards saturation as shown in detail in figure |3] 
At a distance of 1.2/L(m the temperature corrections in- 




200 400 600 800 1000 1200 1400 
Corrugation height a (nm) 



FIG. 3. Effect of the temperature on the Casimir force for 
L=150, 300, 600, 1200 nm (from bottom to top) as a function 
of the trench depth a. 

crease the zero temperature force by ~ 5% between two 
flat plates. Remarkably, this increase becomes ~ 20% if 
the doped Si plate contains deep trenches (a 1.4yLtm). 
For L ~ 600 nm the effect is less pronounced but still 
amounts to an increase from 3% to about 10%. Clearly 
the use of a structured surface increases the thermal 
Casimir force and makes the effect easier to be observed 
at shorter distances. A possible explanation is that the 
nanostructures change the spectral mode density espe- 
cially in the infrared frequency domain, so that thermal 
effects become enhanced, as it has already been pointed 
out in heat transfer phenomena between gratings 16, '17]. 

Next, we turn to assess the validity of PFA. Fig- 
ure |4] shows the deviation from PFA as a function 
of the separation distance L and the corrugation depth 
a. By definition rip{L,a = 0) = 1. For a fixed separa- 
tion distance L, the error made by PFA increases with 
deeper trench depth a. We have 'qp{L^a;T = 300 > 
r]p{L,a;T — OK) for all values of separation distance 
L and trench depth a so that a finite temperature is 
seen to always increase the deviation from PFA. At large 
separation distances, we expect PFA to be valid so that 
limx,_i.oo VfIL, a; T) = 1. In particular, for a fixed corru- 
gation depth a, the functions rjp show a maximum for a 
particular distance L = Lmax as illustrated in figure [5] 
Qualitatively, we can say that this value Lmax ~ ci. More 
precisely we find Lj^^x ~ o, for T = K and L^ax ^ o, 
for T = 300 K. Thus the deviation from PFA increases 
with increasing distances for L < a and decreases with 
increasing distances for L > a. The deviation from PFA 
for T 7^ K shows the same qualitative behavior as for 
T = K. 
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FIG. 4. Deviation from PFA rjp as a function of the separation 
distance L and trench depth a. The lower surface is for T = 
K, the upper one for T = 300 K. 
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FIG. 5. Deviation from PFA rjp at fixed trench depth a — 900 
nm as a function of the separation distance L for T = K 
(bottom) and T = 300 K (top). 



4- Comparison mth experimental data 

We arc now in the position to compare our calculations 
with experimental data from [101 111) . In these experi- 
ments the interaction between a nanostructured silicon 
surface and a gold-coated sphere with a radius of 50 /xm 
was measured. The force was detected at ambient tem- 
perature using a silicon micromechanical resonator onto 
which the gold sphere was attached. As the distance 
between the nanostructured surface and the sphere was 
varied, the change in the resonant frequency of the res- 
onator was recorded. This quantity is proportional to the 
Casimir force gradient dpFs^giL) between the gold sphere 
and the silicon grating. Since the separation distance 
between the sphere and the grating is small compared 
to the radius of the sphere, we can relate this Casimir 
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force gradient to the Casimir pressure Fp^g{L) between a 
plate and the grating as dLFs^g{L) = 2nRFp^g{L). Un- 
der these assumptions, the measured quantity dLFs^g{L) 
when normahzed by its PFA value is identical to rjp, i.e. 
dLF{L)/dLFP^^{L) = tjf{L) where we have omitted 
for simplicity the indices indicating the sphere grating 
geometry. 

In figure |6] we plot r]p{L) for a silicon structure with 
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FIG. 6. Experimental data for the Casimir force gradient 
between a gold sphere and a silicon grating with deep trenches 
(black dots with error bars). The data are normalized to their 
PFA expression. This is compared with our calculations both 
for T = K (blue) and T = 300 K (red). 

deep trenches (sample B in Ref. [10]). This sample has a 
period d — 400 nm, trench depth a = 980 nm and filling 
factor / — 0.478. The experimental data are plotted as 
circles with error bars. Our calculations for T = K and 
T = 300 K are plotted as the blue and red curves respec- 
tively. As first reported in [TD], the measured values of 
r]F {L) between gold and silicon surfaces were found to be 
smaller than predictions using perfect refiectors at zero 
temperature [4] . The inclusion of the material properties 
at zero temperature leads to better agreement [5]. Our 
zero-temperature results presented here are larger than 
those in [5] by about dLF/dLF^^"^ « 0.05. We attribute 
this difference to the replacement of the matrix e by a ~^ 
in the lower-left block of the matrix M in Eq. 



(6b I and 

to solving the differential equations ([7| instead of ( 6a I . 

Let us now discuss the thermal effects. The red line 
in figure [6] plots the calculated results for T = 300 K, at 
separation L from « 150 nm to L « 500 nm. Following 
our previous discussion of Fig. ([5]), since the experiment 
was conducted in a regime where L < a, both rjp and 
dpF / dpF^^^ increase with separation distance. In fig- 
ure [6j the difference between the red line for T — 300 K 
and the blue line for T = K is clearly visible. When 
compared to the experimental data measured at ambient 
temperature, the theory curve calculated for T — 300 K 
gives better agreement than the one at T = K. Un- 
ambiguous demonstration of the thermal contributions 
of the Casimir force in this sample, however, would re- 
quire experimental improvements to further reduce the 
measurement uncertainty. 



So far, the thermal contributions to the Casimir force 
has only been observed at distances larger than l/xm be- 
tween smooth surfaces [3]. At smaller distances, the ther- 
mal effects decrease significantly. As shown in figure [3j 
the thermal contributions to the Casimir force between 
fiat surfaces are expected to be only about 3% at ^ 500 
nm. By replacing one of the surfaces with a grating 
with deep trenches, the thermal contributions increase 
by a factor of 3. Nanostructured surfaces therefore hold 
promise for precise measurements of the thermal Casimir 
force. 

Next, we focus on gratings with shallow trenches that 
were used in |llj . The period was again d = 400 nm, 
but the trench depth was only a = 98 nm. The filling 
factor was approximately / = 0.48. Figure [7] shows the 
measured data points from this experiment together with 
the results of our calculations for this sample both at 
T = 300 K (red curve) and T = K (blue curve). In the 
calculation we take into account the exact trapezoidal 
shape of the corrugation profile via a generalization of 
the formalism presented above [18]. As the range of sep- 
aration distances is the same as in the experiment with 
deep trenches the situation now corresponds to a regime 
where L > a. Therefore r]p and thus OlF/OlF^^^ both 
decrease with increasing distance to reach its asymptotic 
value of 1. Again the theoretical prediction at 300K is 
in good agreement with the measured data. The overall 
temperature effect is less pronounced here than for the 
deep trenches as the trench depth is about a factor of 10 
smaller. 
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FIG. 7. Experimental data for the Casimir force gradient be- 
tween a gold sphere and a silicon grating with shallow trenches 
(black dots with error bars). The data are normalized to their 
PFA expression. This is compared with our calculations both 
for r = K (blue) and T = 300 K (red). 



5. Conclusions 

We have calculated the Casimir interaction between a 
plate and a grating at finite temperature. We find good 
agreement between our calculations for T = 300 K and 
experimental data taken at ambient temperature. Even 



7 



though the experiments are performed at relatively small 
separation distances L < 500 nm, the use of gratings en- 
hances the thermal contributions of the Casimir force. 
Our findings provide an alternative approach to study 
thermal Casimir forces without having to reach separa- 
tion distances of the order of microns. 
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